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INTRODUCTION 



The resummation of certain classes of Feynman diagrams to infinite loop order is a powerful method 
^ \ in quantum field theory. A well known example is the hard thermal loop theory [l| , developed in 
^ the context of thermal field theory, which resums all loop corrections which are of the same order as 
^ tree diagrams, when external momenta are soft. 

r-| ". In recent years, another kind of resummation approach, known as two-particle irreducible (2PI) 
effective action theory, has attracted a lot of attention. In the 2PI formalism, the effective action 
■ is expressed as a functional of the non-perturbative propagator jH, [3|, which is determined through 
^ . a self-consistent stationary equation after the effective action is expanded to a certain order in the 
\^ [ loop or 1/A expansion. This self-consistent equation of motion resums certain classes of diagrams to 
^ ■ infinite order. The classes that are resummed are determined by the set of skeleton diagrams that are 
^ ■ included in the truncated effective action. Numerical studies indicate that the 2PI effective action 
^ . theory is very successful in describing equilibrium thermodynamics, and also the quantum dynamics 
ps) [ of far from equilibrium of quantum fields. The entropy of the quark-gluon plasma obtained in the 2PI 
^ ■ formalism shows very good agreement with lattice data for temperatures above twice the transition 
I> ■ temperature [4] • The poor convergence problem usually encountered in high-temperature resummed 
^ ! perturbation theory with bosonic fields is also solved in the 2PI effective action theory Further- 
^ ] more, it has been shown that non-equilibrium dynamics with subsequent late-time thermalization 
can be well described in the 2PI formalism (see ^] and references therein). The 2PI effective action 
has also been combined with the exact renormalization group to provide efficient non-perturbative 
approximation schemes The shear viscosity in the 0{N) model has been computed using the 
2PI formalism f§|. 

The 2PI effective action theory has its own drawbacks and limitations. When the effective action 
is expanded to only 2-loops, the 2PI effective action is complete. However, when the expansion is 
beyond 2-loops, one must use a higher order effective theory to obtain a complete description (9|. 
Higher order effective theories are defined in terms of self-consistently determined n-point functions 
for n > 2. It has been shown that the 2PI effective action is insufficient to calculate transport 
coefficients for high temperature gauge theories [lOj, but that higher order nPI effective actions can 
be used lllll. 
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The 4PI effective action for scalar field theories is derived in Ref . Il2| using Legendre transforma- 
tions. The method of successive Legendre transforms is used in |9|, ll3|. A new method has been 



developed to calculate the 5-loop 5PI and 6-loop 6PI effective action for scalar field theories |14j . 115 



The 3PI and 4PI effective actions have been used to obtain a set of integral equations from which 



the leading order and next-to-leading order contributions to the viscosity can be calculated |16l . |17 

A lot of effort has been devoted to numerical computations in 2PI theories. For higher order 
nPI theories, progress has been made in understanding the analytic structure of the theory, but 
no numerical calculations have been done. One reason is that the renormalizability of higher order 
nPI theories has not been proven. Another is that numerical calculations for higher order effective 
theories are much more difficult than the corresponding calculations with the 2PI effective theory. 
In this paper, we present some preliminary results from our attack on this problem. We work with 
a scalar field theory in the symmetric phase and use the 4PI effective action truncated at 4-Loops. 
The paper is organized as follows. In sections |TT] and IIIII we define our notation. In section |IV] we 
describe the technique we will use to manipulate divergent integrals. In sections IVl and NH we prove 
the renormalizability of the 4- and 2-point functions, respectively. In section IVlII we show that the 
coupling counter-terms obtained in the previous two sections which renormalize the 4- and 2-point 
functions are the same, to within the approximation order imposed by the truncation of the effective 
action. In section IVIIII we present some results from a numerical calculation in 2-dimensions, and in 
section HX] we give our conclusions. 



II. GENERAL FORMALISM 



We consider the following Lagrangian: 



The classical action is: 



S[4>]=SoM + SM, (2) 
S\nt[v] = /" d^xip^{x) . (3) 



4! _ 

In most equations in this paper, we suppress the arguments that denote the space-time dependence 
of functions. As an example of this notation, the non-interacting part of the classical action is 
written: 

So[^] = \jd'xd'yv{x)[iG-,\x-y)]v{y) ^ V , (4) 



iiU + m^) . 



ip=0 

The effective action is obtained from the Legendre transformation of the connected generating 
functional: 

Z[Ri, R2, R3, R4 = j [d^] Exp[z {S4ip] + Rup + ^i?2¥^' + ^Rs^' + ^R^^^)] , (5) 

i?2, R3, Ri] = ~i LnZ[i?i, i?2; -Rs? Ra\ ? 

5W 6W 6W 6W 
m G, V„ VA^W-R.--R.-- R,- - ft- . 
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Connected and proper green functions are denoted Vj^ and Vj respectively, where the subscript j 
indicates the number of legs. They are defined: 



i5W 



(6) 



The equations that relate the connected and proper vertices are obtained from their definitions using 
the chain rule: 



(7) 



We organize the calculation of the effective action using the method of subsequent Legendre trans- 
forms [oi 13 1 . This method involves starting from an expression for the 2PI effective action and 
exploiting the fact that the source terms and R4 can be combined with the corresponding bare 
vertex by defining a modified interaction vertex. The result is: 

r[0, G, V3, V4 = V + ^Tr In Q-' + '-TtGo'G + r,^,[<P, G, V3, V,], (8) 



G, V,, V, 



X 

4!' 



+ 



A 



'G) + <^2[^,G,V3,V4 



where $2[0; G, V3, V4] represents diagrams with 2 and more loops. 

A key point in the proof of renormalizability is the fact that there are different possible definitions 
of j-point functions for a given j. A general proof of renormalizability must demonstrate that all 
j-point functions are made finite by a set of momentum independent counter-terms which have the 
same structure as the terms in the original Lagrangian. Further, one must show that all j-point 
functions for a given j are made finite by the same counter-term, at a different level of truncation 
of the effective action. Equivalently, since the nPl effective theory coincides with the untruncated 
bare perturbative theory in the limit — )■ 00, the counter-terms of the effective theory and the bare 
theory must also coincide in this limit. 

In this paper we consider only the self-consistent 2- and 4-point functions in the symmetric phase. 
These are obtained by solving simultaneously the equations of motion: 



6T[(P,G,Vs,V,] 



6G 



0. 



5T[<P,G,V,,V,] 



(j,=0,G=G,V3=0,V4=V4 



5Va 



0. 



0=O,G=G,y3=O,V4=V4 



(9) 



From now on we drop the subscript on the 4-point vertex function and write V := V4 and V := V4. 



III. RENORMALIZATION 

The variables m. A, G and V in section [ITl should all carry a subscript B to indicate that they are 
bare quantities. These bare quantities (with subscript B) are related to the renormalized ones by 
the following relations: 

W = Zm% - , 5X = Z^Xb - X , (10) 
Gb = ZG, Vb = Z-^V , 

ZGg^ = Go 1 + 6Gq^ , 6Gq^ = i[bZU + brn^) , bZ = Z-\. 

In order to simplify the notation we have not introduced a subscript R for renormalized quantities 
and suppressed the subscript B everywhere except in the equation (fTOl) . All quantities in section [ITl 
are bare, and all quantities in all other sections are renormalized. 
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We consider the 4-loop 4PI effective action. We introduce a finite set of counter-terms witli 
tlie same structure as tlie terms in tlie original action. Tlie effective action expressed in terms of 
renormalized quantities is: 



iV[G, V]= -^Tr In G'^ - ^TrGo ^G + $i,t[G, V] 
$A = -^Tr5Go ^G 



^(A + 5Aet)G2 + l(A + (5Abb)GV, 



2 4! 48 



Tlie functionals $a and $b are shown in figure [m In all diagrams, bare 4- vertices are represented 
by white circles, counter-terms are circles with crosses in them, and solid dots are the vertex V . The 
notation 5\et and 5\bb indicates the counter-terms associated with the bare vertices on the EIGHT 
and BBALLq diagrams, respectively. The subscripts are introduced only for notational convenience. 
In section IVIII we show that these counter-terms are the same, to within the approximation order of 
the skeleton expansion. 



$R = 




"48 



+5 




EIGHT 




EIGHTc 




BBALLq 



FIG. 1. The effective action at 4-Loops including counter-terms. Tlie open circle is a bare vertex, the circle with a cross denotes 
a counter-term, and the proper vertex is indicated by a solid dot. We label the EIGHT and BBALLq diagrams, which contain 
a bare vertex, and the corresponding counter-terms. 



IV. VERTEX SPLITTING 



We want to remove divergences in the 2- and 4-point functions by adding counter-terms to the 
action. The first step is to manipulate the divergent integrals so they can be written as a sum of 
terms of the form: 

J*-^-* = (finite) + (overall p^independent divergence) (12) 
J*-^-* = (finite) + (overall p-independent divergence) + (overlapping divergences) 

where the notation J*^^^ indicates an a;-loop integral. An overlapping divergence appears in a multi- 
loop diagram when a sub-divergence is present. For the theory to be renormalizable, we must show 
that for all j-point functions, the overlapping divergences can be cancelled by counter-term insertions 
in lower loop integrals, and the overall momentum independent divergences can be absorbed by 
overall counter-terms. 



^ All figures are drawn using jaxodraw llSll . 
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We start by defining some notation. We use the momentum variables: 

Lt = L + K-P, Ls = L + Q-P, Rt = R + K-P, Qt = Q + K - P . (13) 

Since the vertex is symmetric, we can write the four momentum arguments in any order. We write 
all momenta arguments as incoming. We write propagator arguments as subscripts to save space: 
for example, G{L) =: Gl. 

Now we discuss how to divide integrals into the form in equation ( IT2|) . The basic strategy is 
as follows. We write some group of factors in the integrand as a function f{{Pi},L) with L an 
integration variable and L ^ {Pi}- Then we split the function / by writing: 

/({P,}, L) = A^f + /(O, L) with A^f = /({P,}, L) - /(O, L) . (14) 

For large L we have A^f ~ j^f. Using Weinberg's theorem we can show that the integrals we 
need to deal with are at most logarithmically divergent, and therefore a term that contains a factor 
Al/ ~ t;/ is finite when integrated over L. 
Some of the diagrams we will need to consider are shown in figure [2l 




L+K-P 



L+K-P 



P L+K-P K 




Q L Q+K-P 

(b) 




FIG. 2. Some of the loop diagrams that we will study in this paper. 



We illustrate the splitting procedure by looking at the 1-loop diagram in part (a). The integral is: 



iri'P,Q,K,-Qt 



dL V{Q, -Qt, L, -Lt)GLGL,V{-L, L„ -P, K) 



(15) 



We define: fi{{Pi},L) = V{Q, -Qt, L, -Lt)GL and f2{{Pi},L) = Gl,V{-L, Lt, -P, K) and use 
([Til) to write: 



iri~P,Q,K,-Qt 



dL[{ALh + /i(0, L))(Ai/2 + /2(0, L))] . 



(16) 



Expanding the square bracket gives four terms. The three terms that contain at least one Al are 
finite. The remaining term is momentum independent. We have: 

iri-P. K, -Qt) = /r''"°H-P, Q, K, -Qt) + /i(Vo, Vo) , 
where we written the momentum independent divergent term: 



(17) 



h{Vo,Vo[ 



dL V{0, 0, L, -L)GlV{L, -L, 0, 0) . 
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Equation (|T7|) is shown schematically in figure |3l 




FIG. 3. Diagrammatic representation of equation H17[l. The grey blob represents the momentum independent divergent integral 
-^'iCVo, Vb) defined in equation p8|) . 



We will also need to consider 1-loop graphs like the one in part (a) of figure [2] but with one or 
two bare vertices. The procedure is exactly the same. For example, if the upper vertex is bare, we 
define /i = \Gl and repeat the calculation above. Note that A^/i is zero in this case, so there are 
only two terms in the expansion, one of which is finite, and the other momentum independent. We 
write the result for one and two bare vertices: 

/,-^^(-P, Q, K, -Qt) = Q. K, -Qt) + /i(Ao, Vo) , (19) 

I^\-P, g, K, ~Qt) = Q, -Qt) + /i(Ao, Ao) . 

We can follow the same procedure for the u- and s-channels. For example, the integral for the 
s-channel diagram in part (b) of figure |2] is: 

I^i-P, K, -Qt) = j dL V{Q, -P, L, -Ls)GLGLy{-L, L„ K, -Qt, ) . (20) 

We define A = V{Q, -P, L, -L,)Gl and /a = ^^.^(-L, L„ K, -Qt). Splitting gives: 

I^i-P^ K, -Qt) = /i^^"°)(-^, K, -Qt) + hiVo, Vo) . (21) 

One can follow the same procedure for all 1-loop graphs. The finite contribution depends on the 
channel, but the momentum independent divergence is the same for all three channels. We write 
the results as in equations (fT7|) and (fT^ with a subscript t or m or s on the finite part to indicate 
the channel. 



Now we consider the 2-loop graph in part (c) of figure |2l which we call the triangle graph. The 
upper loop is J]^^(— P, R, K, —Rt) (see equation (120|) ). and thus we have: 

Jtri(-P, P, K, -K) = JdR ir{-P. R, -Rt)GRGR,V{-R, Rt, P, -K) . (22) 

Using equation fl2T]) to rewrite 1^^ {—P, R, K, —Rt) as the sum of a finite and divergent piece we 
obtain: 

JtH(-P, P, K, -K) = JdR [Jf ^'^^)(-P, R, K, -Rt) + hiVo, Vo)] ■ (G^Gr, V(-P, Rt, P, -K)) . 

(23) 

The second term in the square bracket contains an overlapping divergence which we will show can 
be cancelled by a counter-term. We extract this term: 

/tri(-P, P, K, -K)\, =: hiVo, Vo)F{K, P) , (24) 



where we have defined: 

F{K, P) = JdL {GrGr,V{-R, Rt, P, -K)) . (25) 

Now we look at the first term in the square bracket in equation (123|) . Defining fi{{Pi},R) = 
li^^^'"\-P,R,K,-Rt) and f2{{Pi},R) = GRGii,V{-R, Rt, P, -K) and using (JUD we get four 
terms. The three terms that contain at least one Ar are finite. The remaining term is a momentum 
independent divergence: 

h4-P, P, K, -ir)|di.erge„ce= / dR /(^^^^^^(O, P, 0, -R)G\R)Vi-R, P, 0, 0) , (26) 

=: Is{Vo,huhfo) , 

where the notation bubjo indicates the finite part of the bubble insertion in the s-channel. Combining 
the results in flMl) and fl26l) we have: 



IU-P,P.K,-K) = finite + h{V^,Vo)F{K,P) + hiVoM'ofo) ■ (27) 
This result is shown schematically in figure |H 



finite + 






FIG. 4. Diagrammatic representation of equation (|27l 



Now we look at the 2- loop diagram in part (d) of figure |2] which we call the double-scoop diagram. 
This diagram is a little harder to handle because of the internal vertex which does not connect with 
any external legs. The upper loop is I^^{—P, K, P, —Rt) (see equation (fT5|) ). and thus we have: 



'dscl 



-P, P, K, -K) = j dR I^i-P, K, R, -Rt)GRGR,Vi-R, Rt, P, -K) . (2^ 



We rewrite the t-loop insertion by doing the same calculation that allowed us to rewrite (IT^ as f lT7|) . 
However, the t-loop insertion in fl28p has external legs that connect to the lower loop, and thus they 
carry momenta that are not external to the diagram as a whole. In order to deal with this vertex, 
we split it in external variables only: 

V{Lt, -L, -Rt, P) = [V{Lt, -L, -Rt, P) - V{L, -L, -R, R)] + V{L, -L, -R, P) . (29) 

This splitting is advantageous because the quantity in square brackets is ~ 1/L when L is large, 
and ~ 1/P when P is large, and therefore acts like either or Ar, depending on which integral 
is under consideration. Using this strategy, it is straightforward to see that equations (fT7|) and (fT8|) 
are modified to become: 

lY^{-P, K, R, -Rt) = J,^^^°(-P, K, R, -Rt) + j dL V{0, 0, L, -L)GlV{L, -L, -R, R) . (30) 



In addition, we split the second factor in ([28]) using ([H]) with fi{{Pi}, R) = GRGRy{-R, Rt, P, -K). 
Combining we have: 

h,,{-P,P,K,-K) (31) 
= jdR [/f^fi"(-P,ir,i?,-i?i) + J dLV{0,0,L,-L)GlV{L,-L,-R,R)j ■ [A^A + /i(0, i?)] . 

Expanding we get four terms. Taking the first term in each square bracket gives: 

termi = j dR I^^^'^i-P, K, R, -Rt) A^/i = finite . (32) 

Taking the last term in each square bracket gives a momentum independent divergence: 

term2 = J dR J dLV{0,0, L, -L)GlV{L, -L, -R, R)GlV{-R, R,0,0) =: l2{Vo,V,Vo) . (33) 

The first cross-term in (13T]) is: 

crossi= [ dR [ dLV{0,0,L,-L)GlV{L,-L,-R,R)ARf^. (34) 



The integral over R is finite, but the integral over L is not. We divide the L integral into a finite and 
divergent piece by splitting the middle vertex. Choosing f2{{Pi},L) = V{L, —L, —R,R) and using 
equation (fT^ we obtain a term containing a which gives a finite result when the L integral is 
done, and a momentum independent divergence. The result can be written: 

crossi = finite + j dR j dL V{0, 0, L, -L)G|V(L, -L, 0, 0) Ar/i , (35) 
= Ii{Vo,Vo) J dR{GRGR,V{-R,Rt,P,-K) - GlV{-R,R,0,0)) , 

= h{Vo, Vo)F{K, P) - h{Vo, Vo)jh{X, Vo) , 

where we expanded A/j/i in the second line, and used (fTSl) and (l25l) in the last line. 
The second cross-term in (13T|) is: 

cross2= / dR /,^^*^"(-P, K, R, -Rt) /i(0, R) , (36) 



jdRj dL[Vi-P,K,L,-Lt)GL,GLVi-L,Lt,R,-Rt)-V{0,,0,L,-L)GlV{-L,L,R, -R)] 
■ GlViR,-R,0,0), 



where we used fl30|) to expand the insertion . The integration over L is finite, but the integration 

over R is not. We split both of the middle vertices to construct a A^ by choosing /3({Pj},i?) = 
V{-L, Lt, R, -Rt) and U{{Pi}, R) = V{-L, L, R, -R). The sum of the two A^ terms is finite. We 
obtain: 

cross2= finite (37) 
+ j dL j dR[V{-P, K, L, -Lt)GL,GLV{0, 0, R, -R) - V{0, 0, L, -L)G|y(0, 0, R, -R)] 
■ GlV{R,-R,0,0), 
= finite + F{-K, -P)h{Vo, Vo) - jh{X, Vo)h{Vo, Vo) , 
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where we used (fTSj) and (125|) in the last hne. The full result for the double-scoop diagram is 
obtained by combining equations (132|) . (133|) . f l35|) and f l37|) . In many equations and figures in this 
paper we do not include separately contributions that correspond to different permutations of the 
external legs, but use numerical factors in brackets to indicate the presence of terms or diagrams 
that correspond to permutations which are not written explicitly. Using this notation we write 
F{K, P) + F{—K, —P) = {2)F{K, P) and the full result for the double-scoop diagram is: 



lU-P, P, K, -K) = finite + (2)F(K, P)Ji(V^o, ^o) - 2-/i(A, ^) + h{V,, V, V^) . 



(38) 



This result is shown schematically in figure [5l The factors (2) in the diagrams indicate that the 
graphs with the grey blob on the lower vertex are also included. 




= finite + (2) 





FIG. 5. Diagrammatic representation of equation (f38 



V. RENORMALIZATION OF THE 4-POINT FUNCTION 

The equation of motion for the 4-point function is obtained from the variational equation 5T /5V = 
0. From this point onward, we suppress the last momentum argument in the 4-point function V, 
since it can be obtained from the conservation of momentum. The effective action in figure [1] gives: 



V{Pa, Pb, Pc)=X + SXbb + Vs{Pa, Pb, Pc) + Vt{Pa, Pb, Pc) + K(Pa, Pb, Pc) 



(39) 



with 



K(Pa, Pb. Pc) = ^J dQ ViPa, Pc, Q)G{Q)GiQ + + P^WiP,, P^, -Q) , 

Vt{Pa,Pb,Pc) = Vs{Pa,Pc,Pb), 
Vu{Pa,Pb,Pc) = Vs{Pa,Pb,Pd), 



(40) 



where P^ = —(Pa + Pb + Pc)- The corresponding diagrams are shown in figure [6l 
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FIG. 6. Self-consistent equation of motion for the proper 4-point vertex. 

Combining permutations of external legs, the diagrams in figure are represented as in figure [3 
The factor (3) indicates that there are 3 channels, only one of which is drawn. 




FIG. 7. The equation of motion in figure [6] with only one channel drawn explicitly. 

From equations f|T7|) and fl2T]) it is clear that all three channels will produce the same divergent 
contribution. Using the result in equation f|T7|) and figure [3], the vertex is renormalized by a 1-loop 
counter-term: 

^aS = -^A(H,H). (41) 

VI. RENORMALIZATION OF THE 2-POINT FUNCTION 

In this section we discuss the renormalization of the 2-point function. In section IVI Al we show 
that the 2-point function can be written in terms of another 4-point vertex which we will call M. In 
section IVlBl we prove that this 4-point function can be renormalized with an appropriate choice of 
the counter-term (5Aet- In section |VTT] we will show this counter-term is the same as the counter-term 
5Abb which was used to renormalize the 4-point function in section |Vl to within the approximation 
order of the theory. 



A. Definition of the Bethe-Salpeter vertex 

We show below that the divergent terms in the 2-point function can be packaged into a 4-point 
vertex. The discussion in this section follows exactly that of reference fl9| and is included only for 
completeness. The self-energy is defined through the Dyson equation and the equation of motion for 
the 2-point function: 

G-^ = Go^-S, S = 2^. (42) 
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Including all momentum variables, the expression for the self-energy is: 

S(j9) = ti6ZP^ - W) + ^(A + 5Aet) j dQ G{Q) 

+^(A + 5Abb) j dQ j dK V{P, Q, K)G{Q)G{K)G{Q + K + P). (43) 

In order to extract the leading asymptotic behaviour we define: 

G = G'as + 5G, £(6-) = Sas + Eo, (44) 

where Gas ~ and Sas ~ contain the leading asymptotic behaviour of the propagator and 
self-energy respectively. From 0421) (and using the bare propagator defined in (|2])) we obtain: 

G = Gas + Gli-tm^ + So) + ■ ■ ■ ^ 6G = GH-im^ + Sq) + ■ ■ ■ . (45) 

We define two new 4-point vertices: a kernel A, and a vertex M obtained from A by solving 
a self-consistent integral equation. Both of these vertices will only be needed with momentum 
arguments (P, —P, K, —K) and we will use the short hand notation A(P, —P, K, —K) = A(P, K) 
and M(P, -P, K, -K) = M{P, K). We define: 

M{P, K) = A(P, + dRA{P, R)G{RfM{R, K) . (47) 
From equations (jH]) and f H6l) we obtain: 

Sas + So=S[Gas + 5G], (48) 

^ ^q{K)= -i6m^ + dP6G{P)A^{P, K) + finite terms , 

= -i5m^ + ^ j dP6G{P)A^{P, 0) + finite terms . 

To obtain the second line we expanded E[Gas + SG] in 6G. In the third line we use the same strategy 
as in section |Vl we write Aas(P, K) = [Aas(P, K) — Aas(P, 0)] + Aas(P, 0), and use Weinberg's theorem 
to show that the square bracket produces a finite integral. 
Next we observe that equation (H71) gives: 

Aas(P, 0) = M(P, 0) - ^ y dRAUP, R)Gl{R)M{R, 0) , (49) 



and substituting ( H9l) into equation (HHi) gives: 

So(/s:)= -i6m^ (50) 
1 f . 1 

Using f H5|) and f HHj) to replace the two square brackets in f l50|) we have: 

Eo(-ft')= —i5rn? + finite terms (51) 



dP[5G{P)]M{P,0) - ^ j dR^j dP5G{P)A^{P, R) GL(P)M(P,0)+ finite terms. 



+ ^y" dP[{-im' + J:o{P))GUP)]M{P,0)-^j rfp[So(P)+«W]GL(P)M(P,0) 
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The key point is that the terms involving Eq on the right side cancel. We are left with: 

So(K) = -i(5m2 - ^(m^ + 5m^) j dPGl^{P)M{P, 0) + finite terms . (52) 

The P-integral is divergent but momentum independent, and can be cancelled by an appropriate 
choice of the counter-term 5m^. Thus we must show that the vertex M can be renormalized. 

B. Renormalization of the vertex M 

In this section we discuss the renormalization of the vertex M defined in equation fl47p . This 
equation is represented diagrammatically in figure [HI 




+ • • • 



FIG. 8. Diagrammatic representation of equation (|47|) . The black box is the vertex M and the light grey box is the vertex A. 

We will define M(^) to mean all terms produced by iterating x times, and m^^^=M^^^-M'^^ means 
terms produced at the x-th interation which were not present at the (a;-l)-th iteration. 

The vertex A is obtained from equations (fTT]) and (H6|) and shown in figure |9l It contains only t- 
and M-channels. The t-channel is drawn, and the factors (2) indicate that the w-channel is included 
although it is not drawn. Diagrams that are not symmetric when inverted top to bottom carry a 
factor (4) to indicate that there are actually four permutations, only one of which is drawn. 




BBALLct Triangle 
Double Scoop 



FIG. 9. The kernel A. 

Consider for the moment the 1-loop contributions to M*^^^ = m^^'' + m^^^ . The t- and ^-channels 
appear in m'^^^ = A. The s-channel is the part of m^^^ which is produced by the 0-loop piece of A: 

= dLA{P, L)G\L)A{L, K) = j dLXG\L)\ + ■■■ . (53) 

Combining these pieces gives the result is shown in figure [101 
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m 



(0) 



m 



(1) 





FIG. 10. 1-loop contributions to the vertex M. 



Thus we find that M^^^ contains 1-loop diagrams which correspond to all three channels, but the t- 
and M-channels are produced at the level of 0- iterations (in the kernel), and the s-channel does not 
appear until the first iteration. It is clear that this structure also holds at the 2-loop level: t- and 
M-channels are contained in the kernel and s-channels are produced by iteration. 

In order to show how renormalization works, we divide the counter-terms into pieces and group 
each piece together with the divergence it cancels. To do this, we introduce a parameter c which 
serves as a counter so that all terms in m^^^ are proportional to c^. The parameters c is introduced 
only for organizational convenience and will be set to one at the end of the day. Using this notation 
we write the counter-term in the Lagrangian: 



5X 



et 



5A^P + AA 



et ) 



AA, 



■et 



(54) 



where the notation AAgf *■* indicates the contribution to the counter-term A A produced by the ith 
iteration. The factor c in the expansion of AAgt ensures that contributions from this counter-term 
are grouped not with terms from the kernel, but with terms that are produced by iteration. This 
grouping allows us to show that all overlapping divergences cancel, and to construct the counter-term 
that cancels the remaining overall divergence. 

As a first step, we substitute equation (15^ into the expression for the vertex A and write the 
result: 



A = A[(5A,t, 5Abb] + AAet =: A^ + AA, 



et 



(55) 



The notation A [5 Ag^^, 5 Abb] means the vertex in figure [9] with the full counter-term 6\ct replaced 
by only the piece SX^^. The counter-terms 6X^f and 5Abb are chosen to make A finite, and this 
finite piece is denoted Aj. These counter-terms can be obtained from figure [HI using the results in 
equations ([IT]), (127]) and ([MD (see figures El H and [S]) . We obtain: 



6X, 



bp(i) 



et 



''^bb 



6X, 



bp (2) 



et 



= -(4)i/i(A,K,) + (2)^/i(K.,^o), 
-^/i(K.,K,), 

= -(2)i [h{Vo, V, Vo) - 2jhiX, Vo)h{Vo, Vo)] - (4)^/3(^0, bub/o) 



(56) 
(57) 
(58) 



where the super-scripts indicate the loop order of each contribution. Note that equations (HTl) and 
f l57)) agree. 



Now we look at the vertex M. We rewrite equation (14 7p with a factor c inserted to keep track of 
the iteration order: 



M(P,K) = Af{P,K) + AX. 



et 



1 

+ 2 = 



dR 



Af{P, R) + AAgt G{RfM{R, K) 



(59) 
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We show below that this equation can be renormahzed using AAet self-consistently determined from 
the equation: 



AAet = -\c j dL [A/(0, L) + c AAet] G\L)M{L, 0) . 



(60) 



Since both equations (l59l) and (l60l) involve only the s-channel, we can suppress all momentum 
arguments without ambiguity if we maintain the order of the factors in every term. Using this 
notation equations (l59l) and f l60l) become: 



M = 

AAet 



A/ + AAet 

1 



1 



A/ + AAet 

2 



[(A; + cAAet)G'M]oo, 



(61) 
(62) 



where the subscript 00 in the second line indicates that the legs on both sides of the diagram carry 
zero momentum. 

Now we substitute AAet = Xli c^^-^cf *^ (see equation ( l54l) ) and collect powers of c. At order we 
obtain: 



c°: AA('*°) = 0, M(o)=m(°)=A 



/ 



At order we have: 



AA, 



(iti)_ 



et 



--[A^G2A^]oo = --/i(A/o,A/o) 



(63) 



(64) 



where in the first line we have used the same notation as in equation fllSp with the proper 4-vertex 
replaced by the vertex Aj. To see that m*-^-' is finite, we split the vertex Aj using the same method 
as in section IIVI We introduce a diagrammatic notation which is shown in figure [TT] 



P 




Af (P, L) Al = [Af (P, L) - A^O, L)l A^O, L) 



FIG. 11. Splitting the kernel A/. The arrow on the grey box indicates that the vertex has been split to the right. The integral 
over the loop connected to the legs on the right side of the box will be finite. 



We define a general procedure for splitting s-channel graphs of the form shown in figure [8] at any 
order in c. In all graphs, start from the right-most Aj vertex and split it towards the left. If the loop 
to the left of this vertex is made finite by the splitting then split the next vertex left, if the loop to 
the left of the original vertex is not finite, then split the next vertex right. Continue moving across 
the diagram towards the left until reaching the left-most vertex. Split this vertex to the right. The 
application of this procedure to the vertex m^^^ (see equation (1M|) ) is shown in figure [T2l 
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2m(i) = 




FIG. 12. Splitting the first iteration of M. Tlie divergent piece cancels. 

Now we look at terms of order c^. The counter-term (from equation (l62|) ) is: 

AAif ) = -l[AA(r^G^A/]oo - \[^fG'm^%o. (65) 

We substitute m*^^) and AA^'*^^ from the previous step. Note that the first diagram in the last line of 
figure [12] drops out because the split vertex on the right end has the form A^ = Af{L, K) — Af{L, 0) 
and gives identically zero when we set the external momentum to zero. In the second diagram we 
re-expand the split vertex and write the full result as: 

AAif ^ = l/i(A;o,A/o)^/i(A,A/o) - ^/2(A/o, A/, A/o) , (66) 

where we have used equation flM|) and the same notation as in equations f|T8|) and fl33|) with the 
proper 4- vertex replaced by the vertex Af. Equation (166|) is shown in figure [T3l 



4AA.". . 2 ycjXZK - >00< 

FIG. 13. Diagrammatic representation of the counter-term AA'''^-* in equation (|65p . 

The vertex m^^-* (from equation ( l6T]) ) is: 

m(^) = UxS'^G^Aj + lA,G^m« + AA^f ^ . (67) 

The graphs obtained by substituting the results for AA^f m*^^-* and AA^f from previous steps 
are shown on the left side of figure [TH The results of the vertex splitting procedure are shown on 
the right. Terms that cancel are marked in pairs with the letters (x), (y), (z), {u) and {v). The 
remaining terms are finite. 
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(ic^^ ^GO>^) 

(;/) (^0 (") 




FIG. 14. Splitting 4mf2' towards the left. 

It is straightforward to see how to apply this procedure at an arbitrary number of iterations. 
Each iteration produces one diagram of the form AjG^A/G^A/G^ ■ ■ '^Z which contains a series of 
overlapping divergences that have the form of the terms on the right side of the first line in figure 
[m The counter-term obtained from (160|) contains exactly the right series of nested insertions to 
remove these overlapping divergences. 



VII. EQUIVALENCE OF THE COUNTER-TERMS 
A. Renormalization of the vertex M to 2-Loops in the skeleton expansion 

In this section we study the renormalization of M to 2-loops in the skeleton expansion. In the next 
section we will show that the counter-terms that make M finite at 2-Loops are the same as those 
that make the proper vertex V finite when it's calculated from the equation of motion to 2-loops in 
the skeleton expansion. 

For convenience we collect all pieces of the counter-term 6Xet which are produced by the first two 
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iterations (see equations (156|) . (158|) . f lM|) and ( 166|) ^ 



+ 



sx 



bp(l) 



et 



5A 



bp (2) 



et 



AA 



(iti) 



et 



AA 



(it 2) 



-(4)i/i(A,K,) + (2)iji(K,,K.), 

-(2)i [/2(K), V, Vo) - 2i/i(A, Vo)/i(K), Vo)] - {^)h,{Vo, bub^o) , 
1 



/i(A/o,A/o) 



J2(A/o,A/,A/o) -2/i(A/o,A/o)^/i(A,A/o) 



(68) 
(69) 
(70) 
(71) 



Recall that 6X^J^^ and ^A^^'' are the 1- and 2-loop counter-terms that renormalize the kernel A, and 



.(2) 



AAi;"^ and AA^""^ are the counter-terms that renormalize M to 1- and 2-iterations. Since the vertex 
Af contains 2-loop terms, we need to separate the 1- and 2-loop pieces of AAit"\ We write: 



vt 

,(it2) 



AA(r^ = AA, 



(iti) 



1 loop 



+AAir) 



2 loop 



AA, 



(iti) 



et 



1 loop 



2 loop 



-2^i(A,A), 



-{2)[AfG\L)A 



(1)1 

/ Joo' 



(72) 
(73) 
(74) 



where Aj*"* and A^p are the 0- and 1-loop parts of the finite kernel. 



The 1-loop counter-term is obtained from equations (168|) and (173|) : 



5Ali) = 5A^f«+AA(^") 



*ct 



= -(4)i/i(A, Vo) + {2)h{Vo, Vo) - ^/i(A, A) 

1 loop z Z Z 



(75) 



Equation (175|) is shown in figure [151 



5X^ 



5A^t^W 



AAir^i 



Hoop 




1 . 



FIG. 15. Diagrammatic representation of the 1-loop counter-term in equation ([75 



Now we look at the 2-loop terms. At 2-loop order in the skeleton expansion we can rewrite the 
2-loop contribution to AAJ^"^ in equation dH]) using Af^ = A ^ V and Af\X,V) Af\v,V) 

where AP{V, V) is the two 1-loop graphs in figure |9] with the bare vertex A replaced by the proper 
vertex V. The result is shown in figure [161 
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(itl), 

et 1 2 loop 




/3(Vb,bubjo) 



FIG. 16. 2-loop contribution from AX^\ 



(itl) 



Similarly we can rewrite the result for AAgf in equation (17T]) by replacing Aj — )■ A^''' = X ^ V. 
The result is: 



AA 



(it 2) 



et 



2 loop 



hiVo, V, Vo) - 2-Ji(A, Vo)h{Vo, Vo) 



(76) 



Combining the 2-loop contributions in equations f TMl) and ( 1761) (see figure [T6l) and comparing with 
equation (1^ we have: 



AA 



(2) 
ct 



AA 



(it 2) 



et 



1 loop 



2 loop 



From equations (169|) and ( !77|) the full 2-loop counter-term is: 



'-6X 



bp (2) 



et 



(77) 



5Al?=5A^rnAA2\ (7J 
= -{3)[^[l2iVo,V,Vo) -2^h{X,Vo)h{Vo,Vo)] + (2)^J3(K), bub/o)} • 



B. Renormalization of the vertex V to 2-Loops in the skeleton expansion 



In this section we show that the vertices V{—P, P, K, —K) and M{—P, P, K, —K) are renormalized 
by the same counter-terms, to within the truncation order of the theory. First we note that the 
counter-term SX^t in figure [15] is equal to the counter-term 5 Abb in equation (HTj) at the 1-loop level, 
since we can replace the bare vertex with the proper one when working at 1-loop. 

At 2-loops these counter-terms will not be the same. This is because the integral equation for the 
vertex V obtained from the n-Loop nPI effective action contains terms up to n — 3 loops, while the 
corresponding equation for Aj contains up to n — 2 loops (these results follow from the fact that 
functionally differentiating a graph with respect to V opens 3 loops, and differentiating with respect 
to G opens 1 loop). Therefore, in order to compare the counter-terms that renormalize these two 
4-point vertices, we must compare the A obtained from the 4-loop 4PI effective action, and the V 
equation of motion obtained from the 5-loop 4PI effective action. 

The 5-loop 4PI effective action in the symmetric theory contains only one additional diagram (T3 |. 
which is shown in part (a) of figure [T71 The corresponding contribution to the V equation of motion 
is shown in part (b) of the figure. 
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1 

128 




-(3)i 



FIG. 17. The 5-loop skeleton diagram in the effective action and the corresponding contribution to the equation of motion. 

From equation ( l38l) we know that there is an overlapping divergence in this 2-loop contribution to 
the equation of motion, and it appears that there is no counter-term to cancel this overlapping diver- 
gence. However, we know that if we expand the equation of motion for the vertex V perturbatively, 
we will generate an infinite set of overlapping divergences and an infinite set of diagrams containing 
the counter-term Abb- To show that the overlapping divergences cancel, we have to rearrange the 
equation of motion in a more convenient way. We use the M equation discussed in the previous 
section for guidance. We rewrite the 1-loop s- and t- u-channels by substituting the equation of 
motion into itself. This is shown in figure [IHl 





— scti 



' XZX - XZX^' XIX^= XZX^' xzxzx 



■■xzxzx 





set2 



FIG. 18. Rearrangement of the 1-loop diagrams in the V equation of motion. 



Combining 1-loop contributions to the equation of motion we have: 



V^^^ = 5AS + seti 



1 loop 



+ set 2 



1 loop 



(79) 



The 1-loop diagrams in seti and set2 are, by construction, the same as the 1-loop diagrams in m^^^ 
(see figure [To]), and therefore we can use SX\2 = ^^^Jt' (see equation (175|l and figure [T5|) to renormalize 
the rearranged V equation of motion at the 1-loop level. 

Now we look at the 2-loop graphs. As before, since we are working at the level of 2-loop skeleton 
diagrams in this section, we can set the bare vertex equal to the proper one in 2-loop diagrams. We 
insert the 1-loop counter-term 5A^t^ obtained from equation ( !75|) (figure [T5|) with Vq = V, which is 
just — 3/2/i(Vo, Vq) or —3/2 times the grey blob in the second diagram on the right side. Adding the 
contributions in figures [TTT b) and [18] and collecting topologies we obtain the graphs in figure [19] 
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-(2)i >OX +1 XZXIX 

FIG. 19. The 2-loop diagrams in the V equation of motion. 

All graphs which contain a grey blob are produced by substituting the 1-loop counter-term into 
a 1-loop topology in figure [THl The graphs with grey blobs in the first two lines come from the 
1-loop t- M-channels (the third graph in the second line of figure [T8l) , and the graphs with grey blobs 
in the third and fourth lines of figure [19] come from the 1-loop s-channel (the second and third 
graphs in the right side of the third line of figure [T8|) . We have artificially divided these counter- 
term contributions into two pieces so that we can see explicitly the pattern of cancellation of the 
overlapping divergences. The coefficients of the vertical and horizontal triangle graphs on the right 
side of figure [19] are obtained directly from figure [18] The coefficient of the vertical double-scoop 
diagram in figure [T^ is obtained from figure [TH] and the t- u-channel parts of figure [TTT b). Figures 
[TTT b) and [TH] contribute — 2^ and 4|i, respectively, and so total factor is |, which we write as (2)^ 
to indicate that the u-channel should be drawn separately. Similarly, the coefficient of the horizontal 
double-scoop is — ;| from figure [TTT b) and \ + j from figure [T8l which gives j. 

The diagrams on the right side of figure [T9]can be divided into finite pieces, momentum independent 
divergences, and overlapping divergences using (I27p and (13 8 p (see figures [1] and [5]) . The overlapping 
divergences cancel and the momentum independent divergences can be cancelled by the counter-term 
5\\^^. Collecting terms we obtain 5\\^^ = 5X),^ (see equation f[78|) ). 

Thus we have shown that if we calculate the vertex M from the 4-Loop 4PI effective action and 
the vertex V from the 5-Loop 4PI effective action, the counter-terms SXet and 6Xhh are the same at 
the 2-Loop level of the skeleton expansion. 

VIII. NUMERICAL CALCULATIONS 

In this section we calculate numerically the 2-point and 4-point functions in 2-dimensions. We 
work in Euclidean space and redefine the variables: 

A = -iXe , 6X = -i5Xe , V = iVe, G = -iGe , S = -iS^ . (80) 
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In this section all variables are Eucledian and we suppress the subscript E. Note that in 2-dimensions 
the coupling constant and 4-point vertex have dimension 2 (and not as in 4-dimensions). 

We will solve the self-consistent equation of motion for the 2- and 4-point functions using a numeri- 
cal lattice method. For comparison, we start by looking at the corresponding perturbative calculation 
in two dimensions. For the 4-point function the diagrams we need are obtained from figure O with 
proper vertices replaced by bare ones. In 2-dimensions there are no ultraviolet divergences in these 
1-loop diagrams, and therefore we can set the counter-term 5 Abb = 0. From equations ( 139|) and ( HOl) 
one can show that the 4-point vertex to order is given by: 



V{Pa, 



Pb, 



■^+87 



dx 



m' 



W? + X(l -x){Pa + Pb 



+ 



+ X{1 - x)iPt + Pc)^ 



+ 0(A^ 



The diagrams we need for the 2-point function are shown in figure 




FIG. 20. 1-loop and 2-loop self-energy diagrams in perturbation theory. 



A simple calculation gives the 1-loop self-energy: 

— ( 7^; In 47r In — 

57r e m 



ii = — (--7£; + ln47r + ln5^), (82) 



where is a mass scale introduced by dimensional regularization. Si^ is logarithmically divergent, 
but the divergence is independent of the external momentum. The 2-loop self-energy is also 
easily obtained as: 

^'"^^^^ = "6(4^^], "^"^io [y + {l- y)x{l - x)]m^ + y{l - y)x{l - x)P^ ' ^^^^ 

Although the 2-loop self-energy S2L is a function of the external momentum, it does not have an 
ultra-violet divergence in 2-dimensions. The inverse propagator is obtained from the self-energy as: 

G~\P)=P^ + m^ (84) 
A2 



dx / dyj — — ^ -— rf^ + 0{X^ 



6(47r)2 7o 7o '^[y^(^i_y^x{l-x)]m^ + y{l-y)x{l-x)P^ 

It can be renormalized by setting 6Z=0 and choosing to absorb the constant divergence in Ei^^ 
(see equation (ITOl) ). 

The structure of the non-perturbative 4PI calculation is similar. The equation of motion for the 
4-point vertex is shown in figure El The loop diagram is finite in two dimensions, and we drop the 
counter-term 5Abb- The equation for the 2-point function can be obtained directly from equations 
f|TT]) and (H2|) . The result is shown in the first line of figure [211 The diagrams can be rearranged 
by substituting the V equation of motion into the vertex on the left side of the sixth diagram. This 
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substitution cancels the 3-loop diagram and produces the result shown in the second line of the 
figure. 



5G 



4 O 4 



+5 



FIG. 21. The self-energy obtained from equations Hll|l and (|42|l . 

The 1-loop self-energy contains one momentum independent divergence from the tadpole diagram 
which can be absorbed by (5m^ . The 2-loop sunset diagram is finite in two dimensions, which means 
that we can choose 6Z = and drop the tadpole and sunset counter-terms in figure [5T] 

We use a N X N two dimensional symmetric lattice with periodic boundary conditions. The lattice 
spacing is a. In Euclidean space, the momentum is discretized: 

2ti N N 
Pi = ni (i = 1,2), Ui = hi,..., — . (85) 

On the lattice, the equation of motion for the 4-point vertex (equations (l39l) and (HOl) ) is transformed 
into: 

V{Pa, Pb, Pc) = -A + Yl [^(^- Q)G{Q)G{Q + Pa + Pc)V{Pb, Pd, -Q) 

+V{Pa, Pt, Q)G{Q)G{Q + Pa + Pb)V{P,, Pd, -Q) 

+V{Pa, Pd, Q)G{Q)G{Q + Pa + Pd)V{Pb, Pc, -Q)] . (86) 
The 2-loop self-energy on the lattice is easily obtained from figure [21] as: 

= ^r4nl E E ^' K)G{Q)G{K)G{Q + K + P) . (87) 

The gap equation in Euclidean space is: 

G-\P) = Go\P) + ^{P). (88) 

Figure [22] shows the 4-point vertex and the self-energy as functions of the coupling strength A, 
calculated from 4PI effective action and perturbation theory, where we choose the external momenta 
to be vanishing. The lattice spacing a = 27r/(A^m) is adopted and two values of lattice number 
= 16 and N = 8 are used. When the coupling strength is small, for example when A is less than 
2m^, the results obtained from 4PI calculations are consistent with the perturbative ones, for both 
the vertex and the self-energy. However, when the coupling constant becomes large, the perturbative 
results are quite problematic, as shown in the left panel of figure [221 where the 4-point vertex changes 
sign at large A. This behaviour does not occur in the non-perturbative result. 

In figures [23] and [Ml respectively, we show the dependence of the 4-point vertex V{Pa = P,Pb = 
0, Pc = 0) and the self-energy on the external momentum. We look at several values of the coupling 
constant. Calculations are performed both in the 4PI formalism and in the perturbation theory. 
We observe again that the 4PI calculations are consistent with the perturbation ones when the 
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FIG. 22. (color online). Left panel: Comparison of the 4-point vertex with external momenta vanishing as a function of the 
coupling strength A, obtained from the 4PI numerical study and perturbative calculations. For the 4PI calculation, we use 
lattice spacing a — 2Ti/(Nm) and two values of the lattice number A'^ = 16 and N = 8. Right panel; Self-energy with external 
momentum zero as a function of the coupling strength. 




I 1 I 1 I 1 I 1 I 

2 4 6 8 

P (m) 



FIG. 23. (color online). Dependence of the 4-point vertex V{pa — p,pb ~ 0,Pc ~ 0) on the external momentum. Results are 
shown for perturbation theory and using the 4PI effective action formalism. We use lattice spacing a — 2n /{Nm) and A*" = 16 
for the 4PI calculations. 



coupling is small. However, the discrepancies between the two calculations become more and more 
pronounced as the interaction strength increases. 

As a final comment we note that the full momentum dependence of the 4-point vertex can be 
obtained in the 4PI effective action formalism from equation (l86l) . In figure [251 we give a contour 
plot of the 4-point vertex which depicts the dependence of the vertex on two momentum components. 
One can see that the vertex has a minimum at the origin of the coordinates, and the gradient varies 
with a change of direction. 
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FIG. 24. (color online). Dependence of the self-energy on the external momentum with A — 0.5m^ (left panel) and A = 5m^ 
(right panel). Results are shown for perturbation theory and the 4PI effective action formalism. We choose lattice spacing 
a = 2-K/{Nm) and A'' = 16 for the 4PI calculations. 




P,,(m) 

FIG. 25. (color online). Contour plot of the 4-point vertex — V(pai,Pbi) with other momenta vanishing. We use A = 5m^, 
a = 27r/(iVm) and N = 16. 

IX. SUMMARY AND OUTLOOK 

In this paper we have demonstrated the renormahzabihty of the self-consistent 2- and 4-point 
functions obtained from the 4-Loop 4PI effective action. We have shown that the counter-terms that 
renormahze these functions are the same, up to the level of the truncation of the effective action. The 
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equation of motion for the proper vertex V is symmetric in the three Mandelstam channels, whereas 
the vertex A contains only the t- and u-channels. The divergences in the 2-point function can be 
packaged into another vertex M which resumms A in the s-channel. The vertex M calculated from 
the 4-loop 4PI effective action has the same structure as the vertex V calculated from the 5-Loop 
4PI effective action, and the two vertices are renormalized using the same counter-term 6X on both 
the EIGHT and BBALLq diagrams in the original effective action. 

As an indication of the usefulness of the 4PI effective theory, we have solved the integral equations 
in 2-dimensions using numerical lattice calculations. We have compared our results with those 
obtained in perturbation theory when the coupling strength is small and perturbative calculations 
are reasonable. The two calculations are in good agreement in the perturbative region. When the 
coupling is large, perturbation theory is not applicable, and the results of the perturbative calculation 
are not sensible. The 4PI effective action theory is a non-perturbative method which can be used 
consistently when the coupling is large. The full-momentum dependence of the 4-point vertex is 
easily obtained in the 4PI formalism. 

We have only performed numerical calculations in two dimensions, and it would be interesting to 
extend the calculation in this paper to four dimensions. It would also be of interest to calculate a 
physical quantity like the pressure or shear viscosity using the 4PI effective theory. Work in these 
directions is in progress. 
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